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Abstract Resonant absorption and mode conversion are both extensively studied 

mechanisms for wave "absorption" in solar magnetohydrodynamics (MHD). But are 

they really distinct? We re-examine a well-known simple resonant absorption model 

in a cold MHD plasma that places the resonance inside an evanescent region. The 

^— s normal mode solutions display the standard singular resonant features. However, 

rsj these same normal modes may be used to construct a ray bundle which very clearly 

i undergoes mode conversion to an Alfven wave with no singularities. We therefore 

conclude that resonant absorption and mode conversion are in fact the same thing, at 

least for this model problem. The prime distinguishing characteristic that determines 

which of the two descriptions is most natural in a given circumstance is whether the 

converted wave can provide a net escape of energy from the conversion/absorption 

region of physical space. If it cannot, it is forced to run away in wavenumber space 

instead, thereby generating the arbitrarily small scales in situ that we recognize as 

fundamental to resonant absorption and phase mixing. On the other hand, if the 

converted wave takes net energy way, singularities do not develop, though phase 

<-^h mixing may still develop with distance as the wave recedes. 
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C^ !• Introduction 

I The conversion of MHD waves from one type to another has attracted much attention 

^. over many years. Recently, Schunker and Cally (2006) have adapted the dispersion- 

relation-based generalized ray theory formalism of Tracy, Kaufman, and Brizard 
(2003) to the problem of conversion between fast and slow waves, finding that the 
"fast-to-slow transmission" occurring near the Alfven/acoustic equipartition level 
depends crucially on the attack angle made by the incident wave vector to the 
magnetic field. Comparison with exact solutions in an isothermal atmosphere shows 
that this performs very well indeed at small to moderate attack angles (Hansen and 
Cally, 2009). Fast-to- Alfven conversion has not yet been satisfactorily described in 
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this manner, though a wave-mechanical approach shows that it can be very significant 
in certain attack directions (Cally and Goossens, 2008). 

The Generalized Ray Theory (GRT) of Tracy, Kaufman, and Brizard (2003) and 
Schunker and Cally (2006) reduces the problem of mode conversion/transmission to 
a local analysis of a saddle point in an appropriate phase space. This can be done in 
a quite general way, allowing the development of a simple recipe for calculating the 
transmission and conversion coefficients. 

Parallel and seemingly disconnected theoretical developments in MHD wave the- 
ory have focussed on the well known phenomenon of resonant absorption, in which 
a wave is "absorbed" on a certain resonant surface (see Goedbloed and Poedts, 
2004 for extensive discussion). Applications to solar physics are numerous, ranging 
from theories of p-mode absorption by sunspots (Hollweg, 1988), to mechanisms for 
thermalization of Alfven waves in coronal loops (Ionson, 1978) and the damping of 
TRACE loop oscillations (Ruderman and Roberts, 2002). The physical mechanism 
by which resonant absorption operates through phase mixing to absorb surface waves 
is brought out clearly by Cally and Sedlacek (1992). 

But are mode conversion and resonant absorption always distinct processes? In 
this paper we look at perhaps the simplest resonant absorption model in which a 
fast wave approaches and is absorbed by an Alfven resonance in a cold MHD plasma 
(plasma-/? = 0). Now, an "Alfven resonance" can also be interpreted as just an 
Alfven wave running along a particular field line (or flux surface), and so can also 
be thought of as a result of mode conversion. Our aim here is therefore to better 
understand the relationship between the two theories. In this way, we hope to build 
a bridge between mode conversion and resonant absorption theory, or at least to 
establish a bridgehead. 



2. Modal Model and Governing Equations 

Consider a cold (sound speed zero) magnetoatmosphere with uniform magnetic field 
B in the z-direction (cartesian coordinates (x, y, z)), but where density p varies with x 
only. Assume an ideal magnetohydrodynamic linear wave with plasma displacement 
vector of the form £(x,y,z,t) = (£, x (x),£, y (x),Q)exp[i(k v y + k z z — ut)]\ the only 
restoring force is perpendicular to B, so there can be no £ z . The governing differential 
equations are then 

92 , ^ ,A c _ iU ^k 

" v dx' 



(1) 



ky k z I £ y — ifc. 



v dx 1 



where a{x) is the Alfven speed, showing that £ x and £ y decouple to independent fast 
(£ x ) and Alfven (£ y ) waves if k y = 0. The associated WKB dispersion relation is 

v =(^-ki-ki-e)j(^-e)j=^ ( 2 ) 

where k x is the WKB wavenumber associated with x. The two factors represent the 
fast and Alfven modes respectively. 
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Equations (1) may be easily combined into a single second order ordinary differ- 
ential equation for the dilatation x — V- £ = d^ x /dx + i k y £ y (see also Equation (1) 
of Hollweg, 1990): 

2 d fld X \ (3cos 2 6- "- 2a 



^{p£) + - — n — x=0 ' (3) 

where (3 = (3(x) = co 2 /(a 2 k 2 ) — 1 = Q. 2 sec 2 — 1 > — 1, k y = k ± sin0, k z = k ± cosd, 
and f2 = cj/ak ± . Here, k ± = \Zk 2 + k 2 = e _1 is the component of the wavevector 
perpendicular to the direction of stratification x. The corresponding components of 
displacement are £ x = —k~ 2 [3~ 1 x l = — e 2 /? -1 sec 2 6 \' and £ y = —\f3~ 1 k y k~ 2 x = 
— \efi~ 1 sin 6 sec 2 9 %■ The associated wave energy (Poynting) flux is 

„ „ fe 2 sec 2 6 T . , „. e . , 9 ,, ., cos^_ l9 \ , A . 

F = F l lm( X 'x ), ^sm6»sec 2 6»|x| 2 ,-^|$| 2 J , (4) 

where Fq = ujB 2 //i = lj pa 2 , B = |B| is the background magnetic field strength and 
/i is the magnetic permeability. Note that _fo has the dimensions of a flux per unit 
length, or in other words an energy density per unit time. By energy conservation 
(or by reference to the Wronskian of the solutions of Equation (3)), the x-component 
F x is constant, except for a discontinuity at /3 = which is the signature of resonant 
absorption. The corresponding wave energy density is 

E=y>[(co 2 + a 2 k 2 M 2 + a 2 \ X \ 2 ] . (5) 

We shall assume that /? decreases monotonically with x. Assume a rightward 
propagating fast magnetoacoustic wave classically reflects at X = 0, i.e., /3(0) = 
j3 Q = k 2 /k 2 = tan 2 8, and that there is an Alfven resonance at x = X > 0, i.e., 
P(X) = 0. An incident fast wave in Region I (x < 0) reflects in Region II (x « 0), 
but partially tunnels through Region III (0 < x < X) to be resonantly absorbed in 
Region IV (x « X). There remains only an evanescent tail in Region V (x > X). If 
k y = (0 = 0), Region III collapses to zero width and Regions II-IV coalesce. 

With these points in mind, and using the series expansion (39) from Appendix A, 
the total z-component of wave energy flux may be shown to be 

F x = tt hF ta,n 2 6 \ X (X)\ 2 U(X-x), (6) 

where h = —1/(3'(X) > and IA is the unit step function, x i g everywhere continuous 
so there is no need to distinguish between X + and X~ here. 

Equation (5) may be used to define a net energy transport velocity V = F/E. 
Employing the Frobenius series again it is easily confirmed that V = (0((x — 
X) 2 )IA(X — x), 0(x — X), a + 0(x — X)) near the resonance X (where |£| — > oo and 
X is finite). Note that it takes an infinite time to reach the resonance at speed V x . 



3. Mode Conversion 

Mode conversion is the process whereby one clearly identifiable wave type changes 
into another type in some (normally) compact conversion region where the distinc- 
tion between the modes breaks down. However, away from this region, the incident, 
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converted and transmitted waves are easily identified, and conversion/transmission 
coefficients assigned. The precise mathematical description of this process and ap- 
proximate calculation of conversion coefficients may make certain assumptions to 
render the problem tractable, though the wider concept of mode conversion does not 
depend on these. 

One such description is the so-called modular theory developed over the last two 
decades by many authors, summarized by Tracy, Kaufman, and Brizard (2003), and 
adapted to MHD by Cally (2006) . It fundamentally assumes the wavelengths involved 
are large enough to fit several into the interaction region, a point which limits its 
applicability in the current context, and to which we shall return later in this section. 

Modular theory is based on an action principle with action A = J Cchtdt and 
Lagrangian density C = £ D£, where the superscript "ff" denotes the Hermitian 
transpose. For the problem at hand here, the dispersion tensor is simply 



(k x + k z 



^x^y 



h ,.)• h ,, 



(k + k z 



(7) 



Once again we assume a 2 = B 2 / up increases monotonically with x. 

Fixing k y and k z , contours of the dispersion function T>(x,k x ) = det D may be 
plotted in the (x,k x ) plane, displaying a saddle at the point (x*,0), where uj 2 /a 2 = 
fc 2 + \k 2 i.e., P+ = | tan 2 6 — 1. As a concrete example (see Figure 1), we let 
Q 2 = uj 2 /a 2 k 2 = e~ x , for which X = — \n{k 2 /k 2 ) = In sec 2 9, and the saddle point is 
located at a;* = — ln((fc 2 + |fc 2 )/fc 2 ), k x = 0. 
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Figure 1. Dispersion diagram for the case k y = 3, k z = 9 (see Section 5) with u) 2 ja 2 kj_ = e~ x , 
for which X = ln(10/9) and xq = ln(20/19). The heavy full curves correspond to T> = 
and represent respectively the fast (left) and Alfven (right) modes. The dashed curves are 
T> = T>(Xi,, 0), which pass through the saddle point. A few other values of T> are shown dotted. 
The arrows indicate the direction of energy propagation in the fast mode. The shaded region 
between the resonance and the principal saddle axes has area 0.31. 
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Now, we may decompose D = SAS , where A = diag[(a; 2 /a 2 ) — fc 2 , (cj 2 /a 2 ) — fc 2 
contains the two eigenvalues of D, k = k|, and 



L / Ky Kj 



""n \ "*X l^y 



(8) 



where k p = \/k x + k 2 has been introduced for convenience. The two columns of S 
are the normalized eigenvectors of D. 

The mode conversion formalism of Tracy, Kaufman, and Brizard (2003) requires a 
different form though, in which the diagonal elements of the 2x2 dispersion matrix 
represent the asymptotic behaviour of the two mode types far from the interaction 
region. To that end, we introduce 



D 



1 h.2 _ h.2 _ 1 J, /h.2 i OI.2 1 1,2 

2 K P K z 2 Kx V x+ V 2 K V 

1 b 2 " 2 — 1 b 2 — b 2 -U 1 b lb 2 4- 9i-2 

2 K y ~a? 2 K p K z + 2 Kx V x + V. 



D, 



2 



71 ^RAR T , 



(9) 



where 



Y ^ P k x y'fcp + k y J k p + k x yj/k p + k y 

^^ {^k 2 p + k X ^kjTWy ^k 2 p ~k X ^kj+Wy 



(10) 



Both S and R are orthonormal matrices. Here, the loci D a = and Dt, = 
correspond to the level curves of V that pass through the saddle point, i.e., the 
dashed curves in Figure 1, with D a the incoming branch and _D{, outgoing. The off- 
diagonal terms 77 = ky/2 are the coupling coefficients. As anticipated, there is no 
coupling when k y = 0. 

The point of all this is that we have now introduced a more suitable state vector 
U = Q£ where Q = SR in place of £. With it, the Lagrangian 

C = £ H D£ = U H DU (11) 

remains unchanged in form though with the new dispersion tensor D. The compo- 
nents of U represent the amplitudes of the two decoupled (orthogonal) modes with 
eigenvalues D a and D/, in the conversion region. That a transformation leading to 
D is possible follows from Sylvester's Theorem (see Schunker and Cally, 2006 for 
details) and does not actually require explicit construction. 

The classic "transmission coefficient" for energy "transmitting" from the fast to 
the Alfven mode T = t 2 may then be calculated, where 

r = exp(-7r|r/| 2 )* , (12) 

with 77 = j]/!/!?! 1 / 2 and 

B -I D CnV - dD ° dD * 9D a dIh_ 2 2 , kj + kl 

B-{D a ,D b }- ^ Qx ^ Qkx - 2k ± il ^ fT ^ (13) 
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is the Poisson bracket of the two uncoupled dispersion functions. The subscripted 
star indicates that T should be evaluated at the star point x = a;*, k x = 0. With 
the exponential profile f2 2 = e~ x ' h , and introducing the dimensionless parameter 
a = (k ± h) 2 / 3 sin 2 0, it is then easily shown that 



exp 



ttV? ^ 3/2 



1 + cos 2 6 



(14) 



Unfortunately, this result cannot be taken at face value. As stated earlier, modular 
theory relies on the solution being substantially wavelike (WKB) in the interac- 
tion region. The analysis of Cally (2005) makes clear how this operates in terms 
of stationary phase integrals; two neighbouring modes propagating in phase may 
efficiently exchange energy. A practical measure of whether modular theory applies 
is the (dimensionless) area A s in x-k x space of the region bounded by D a , Db, 
and the resonance (see the shaded region in Figure 1). It requires that A s 3> 2w 
(Brizard, Kaufman and Tracy, 2008), which fails spectacularly in the case of Figure 
1 (A s = 0.31). Nevertheless, comparison with numerical solutions in Section 4 will 
reveal that T always provides an upper bound on the wave-energy transmission to 
an Alfven wave. 

Lest the reader think that different or improved expressions for T may be obtained 
by alternate choice of variables (e.g., arc length along a ray and the associated 
wavenumber), it should be born in mind that the independent variables x and k x in 
Equation (13) must be canonical coordinates in the Hamiltonian sense, and that the 
Poisson bracket (and hence T) is invariant under arbitrary canonical transformations. 
The area A s is also invariant under canonical transformation, by Liouville's Theorem. 



4. Numerical Calculation of the Reflection and Absorption Coefficients 

As befits a local analysis, we choose to approximate the arbitrary Alfven profile by 
a suitable simple analytic form which is expected to suffice in general since mode 
conversion is indeed typically restricted to a compact region. The most convenient 
choice is to assume the density scale length h is constant (its value at the resonance), 
whence Q. 2 = e~ x l h and /3 = e~ < - x ~ x ^ h — 1. Hollweg's 1990 choice of a linear profile 
is problematic, as it unphysically results in Q, 2 < for sufficiently large x unless 
truncated by some ad hoc and troublesome boundary. Unfortunately, this simple 
profile fails to produce an equation with known exact solution (as does the linear 
profile). 

Introducing £ = k ± h~ 1 / 3 (x — X) as the new dimensionless independent vari- 
able, and the dimensionless parameters a = (fc x /i) 2 / 3 sin 2 6 and v = (k ± h)~ 2 / z = 
a~ 1 sm 2 6, Equation (3) takes the form 

1 dX ' ' cos 2 6»H =— ) x = (15) 



dcwodc; v «(c) 

where u = (1 — e~ v ^jv. The even simpler linear profile equation 
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Figure 2. Top panel: Absorption coefficient A as a function of 

a = e- 2 / 3 ft 2 / 3 sin 2 e = (k ± h) 2 / 3 sin 2 6 for U 2 = e~ x / h and 6 = 5°, 10°, ..., 85° (top 
to bottom). The heavy dashed curve is the RAT index (i.e., 9 = 0°). The heavy vertical line 
is at <t = erg. Bottom panel: Phase retardation Aif> between incident and reflected waves for 
the same cases, with corresponding colours, i.e., 9 = 0°, 5°, . . . , 85° (bottom to top). 



is recovered in the limit v — > with a held fixed (and therefore sin# — > 0). Un- 
fortunately, neither equation admits exact solution in terms of special functions 
(neither Mathematica 7 nor Maple 12 can solve them, nor do they appear in the 
compendious tome of Polyanin and Zaitsev, 2003). We shall term Equation (16) the 
Reflection-Absorption/Transmission (RAT) equation, and Equation (15) the Expo- 
nential Reflection- Absorption/Transmission (ERAT) equation. Notwithstanding the 
above protestations about the realism of the linear profile, the RAT equation is useful 
here as the 6 — > 0, k ± — > oo, a finite limit of the ERAT equation. 

Figure 2 summarizes the results of a numerical survey of the solutions of the RAT 
and ERAT equations with < a < 3 and 9 = 0° (RAT), 5°, . . . , 85° (ERAT; see 
the appendices for details). Both the total absorption coefficient A and the phase 
retardation of the reflected wave A(f> are plotted. They are also tabulated in Tables 
1 and 2. Strikingly, absorption never exceeds 50%, and vanishes as a — > and 
a — > oo. These absorption curves are similar to those obtained by Hollweg (1990) for 
a two-piece linear-constant Alfven profile and a nearby "vacuum-boundary". 

The modular theory transmission T is plotted against a and 9 in Figure 3 and 
compared with the numerically derived values of !A~ It is obviously totally wrong 
at small a, though broadly correct (in absolute value and shape) at larger a, both 



convRev.tex; 13/07/2010; 0:17; p. 7 



P.S. Cally and J. Andrics 



in terms of a and 6 dependence. Significantly though, it is always an upper bound, 
T > A. 




Figure 3. Left: Transmission coefficient T as given by Equation (14) for 9 = 0° and < a < 3 
(full curve) compared to the numerically calculated absorption coefficient ■%.(&, 0) (dashed 
curve). Right: T and S\.(a, 0) for a = 1.5 and 0° < 9 < 90°. 

The limited applicability of Equation (14) is no refutation of the mode conversion 
concept per se, merely a reflection of the limitations of the basic modular theory for 
this difficult problem involving a resonance. Nor does it render infeasible calculation 
based on conversion theory, since we have in fact determined exact (numerical) values 
for the absorption coefficient -#(<x, 9) for the exponential profile, albeit in graphical 
and tabular form. Interpolation of Table 1 is as good as having an exact formula. 

4.1. Phase 

Turning now to the reflected wave (the "converted" wave in the terminology of Tracy, 
Kaufman, and Brizard, 2003), the "wave conversion coefficient" 



V / 2~7T : 



rjT(— i \r)\' 



(17) 



is such that the (wave-energy) conversion coefficient C — \c\ 2 is exactly 1 — T. 
However, c, which is complex, also contains information about the phase change 
from the incident to the converted wave: the phase retardation due to mode con- 
version is A(f> c = — arg c. This is plotted in Figure 4, and disappointingly bears 
little resemblance to A</> as depicted in Figure 2. This is probably because of the 
confounding effect of the additional reflection component to that accounted for by 
T. 

Phase change is of relevance to seismology, especially as regards the determination 
of wave travel times (Cally, 2009a; Cally, 2009b). 



5. Superposition and Ray Bundles 

We postulate that "trappedness" of the resonance is fundamental to our interpre- 
tation of the physical process in operation: is it mode conversion or is it resonant 
absorption? Our current model does not at first sight allow us to explore the con- 
sequences of opening the resonance, because the resonant field line is being excited 



convRev.tex; 13/07/2010; 0:17; p. 8 



Resonant Absorption as Mode Conversion? 















^^^3^ 


gg<5g^^ 


105° 










10P . 














vy>>w< 
















/ 




§a\V 


100° 














v\\\\ 








M0y 


w 






m 


95° 






V y 








u* 




' y 










\\\ 




J/ 












V\ 
















\ 

















0.0 



0.5 



1.0 



1.5 



2.(1 



2.5 



3.0 



Figure 4. Phase retardation due to mode conversion as predicted by Equation (17). The line 
styles correspond to those of Figure 2, i.e., 8 = 0°, 5°, 10°, . . . , 85° (bottom to top on a < 1.8), 
with the heavy dashed curve corresponding to the RAT case 9 = 0°. 



uniformly along its entire length. But this z-invariance can be avoided by inverting 
the z-Fourier transform to construct a discrete ray bundle that meets the resonance 
only in a limited region. 

Before synthesizing a ray in this way, we should first define an absolute x-scale. 
Previously, the x axis has been translated so that the reflection point is at x = 
for each mode. But this reflection point position depends on k z . For this calculation, 
the zero of x shall be arbitrarily placed at the reflection point for a specified fc z o 
(the central wavenumber of a distribution); for other k z it is at x = Aj^, where 
Afe 2 = In (sin 2 9/ sin 2 9 ) for our standard exponential density profile with h = 1, 
and 9 = arctan(k y /k z ). Similarly, the resonance is shifted to X(k z ) = lnsec 2 6* + 
Afc 2 . In this section we denote the normal mode solution by x(x, fc z ), equivalent to 
x(x — Afe^ , k z ) of the previous sections, and reserve \ f° r the physical space solution, 
though the t and y Fourier transforms are still in operation. Indeed, we still fix a 
single k y and U), but now integrate over k z . With this in mind, the physical space 
dilatation x is 



X{x,z) 



1 

2tt 



e lk > z A(k z )x{x,k z )dk z 



(18) 



In Equation (18), A(k z ) is a compact amplitude which shall be centred on a particular 
k z = k Z Q. If it is a delta- function, we recover a single mode, but by broadening the 
profile in Fourier space we sharpen it in z-space. A ray bundle may therefore be 
constructed that impacts the (spread) resonance layer and may mode convert to 
an Alfven wave that can escape from the interaction region by propagating in the 
z-direction. 

To this end, a gaussian amplitude distribution 



A(k z 



2-7T H exp 



\{k z -k z v) 2 H 2 



i {k z - k z0 ) z 



(19) 



is adopted that has inverse Fourier Transform exp[— (z — zq) 2 /2H 2 + i z k z o]. This is 
centred at k z g in Fourier space and Zq in physical space, and has widths H~ l and H 
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respectively. The incoming component of x is normalized to 1 at the left boundary 

X = XL- 

The inverse transform is obtained numerically using the fast discrete method of 
Bailey and Swarztrauber (1994) based on the Fractional Fourier Transform (FRFT). 
A large number m of points (typically a power of two or other product of small 
primes suitable for FFT evaluation) is used in both Fourier and physical space lying 
uniformly in [k z a — K, fc z o + K) and [— Z, Z) respectively, where K is large enough 
that A(k z ) is effectively zero outside the interval and Z similarly encompasses the ray 
bundle in physical space. The simple piecewise-constant integration rule is implicit in 
the inversion as applied here. Each of x, dx/dx, dx/dz, and d 2 x/dxdz is evaluated 
in Fourier space and then mapped to physical space by Fourier inversion in order to 
construct the required displacements and fluxes. 

Exact integration of the ray equations in the exponential atmosphere f2 2 = e~ x 
yields the ray path x = 21ncos[|zsec#] for a ray passing through the origin. In 
particular, the "skip distance" in the z-direction for this ray is Az = 27rcos#, which 
allows us to estimate a priori the z-range required in the numerical Fourier inversion. 

The following parameters are adopted by way of example: m = 2 8 = 256, k y = 3, 
k z0 = 9, K = 8, H = 0.5, .i L = -5, z = — k, and Z = 5. For this case, 0.2 < a < 1.94 
and 10.0° < 6 < 71.6°, with a = 0.448, 6 = 18.43°, X = 0.105, and A = 0.464 at 
k z = k z0 , suggesting strong absorption. Figure 1 depicts the phase space structure 
for the central wavenumber. For the bundle, the resonance point x = X(k z ) is 
now widely spread because of the distribution in k z space; for |A(fc z )| greater than 
half-maximum it spans the interval —0.35 < x < 0.7. 

Figure 5 displays the result of this calculation as contour plots of Rex and \x\ 
over the x-z plane. A clearly defined ray bundle is apparent which fits very well 
with the analytic ray path (in red). It enters around z = —it on the left, reflects 
near x = 0, and leaves the domain at positive z on the left. This is as we might 
expect. But note that the intensity of the bundle is less by some 26% on exit than 
on entry, indicating the energy must have gone elsewhere. We should rather look at 
F z = F Re{i$-de/dz}. 

Since x is a C 1 function Fourier inversion results in x being C 2 . On the other hand, 
the singularities in £ make direct numerical inversion problematic. It is preferable 
instead to construct the displacement from x by solution of ordinary differential 
equations in physical-space: 

d 2 £ LJ 2 
az 1 a 1 

where V p is the x-y part of the gradient operator. Equation (20) is easily derived 
from Equations (1), and is solved subject to £ — » as z — > — oo, yielding C 3 solutions. 
The interpretation of Equation (20) is clear: £ propagates Alfvenically along the field 
lines with — V P x acting as a source term. The exact solution is 

1 f z 
£(x,z)=- / 'VpxiXyz') sin k(z' — z) dz' 

X (21) 



^-J_J P xM 



e ifc(z'-z) _ e -ife(z'-z) 



dz' . 



where k(x) = co/a = f2fc xo , and k ±0 = \/k 2 + k 2 . Resonant transfer of energy from 
X to £ clearly occurs where V P x has a substantial Fourier component resonant with 
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Figure 5. Contours of Rex (left) and \\\ (right) constructed by inverse Fourier Transform 
plotted in the x-z plane for the example case discussed in the text of Section 5. The grid used 
is 561 X 256. The exact ray path of the central k z is over-plotted in red, and the green line 
indicates the position of the resonance for this central value. The ray enters the domain from 
the lower left boundary and leaves, with reduced intensity, via the upper left boundary. Dark 
purple in the right panel corresponds to and maximum values to white shading. The contour 
levels are evenly spaced. 



k. In practice, the integrals are performed on the grid using a modified trapezoidal 
rule 1 and the lower boundary condition is applied at z = — Z on the basis that there 
is no significant signal below that point. 

Equation (21) may be interpreted more generally. The sine, or the combination 
of exponentials, may be thought of as eigenfunctions of the Alfven operator along 
the field lines. The displacement £ is therefore constructed by projection of the 
compressional source term onto these eigenfunctions. We believe this result will carry 
over to more complex examples than the simple case considered here. 

Though x is a poor marker of an Alfven wave, £ y is ideal. It is displayed in the 
top left panel of Figure 6, making very clear the generation and propagation of an 
Alfven wave in the resonant slab. Phase mixing is also apparent as the Alfven wave 
propagates upwards, caused by differing propagation speeds on neighbouring field 
lines. Now, the x and z components of wave energy flux are F x = F R- e (ix* 6x) an d 
F z = F Re(£, x d£,*/dz + £ y d£*/dz) = F Ztfast + F ZtAV , where F zAU is that part of 
F z associated with £ y . The top right panel of Figure 6 again convincingly shows the 
presence of an Alfven wave taking energy upwards and away from the interaction by 



1 The trigonometric form is used, with sinfc(z' — z) first expanded by multiple angle formula. To 
ameliorate the effects of imperfect resolution of cosfcz' and sinfcz' on the grid z\, Z2, ■ • • , z m , 
we apply a piecewise-linear approximation for V P x an d perform the resulting trigonometric- 
by-linear integrals exactly. Higher order schemes are unnecessary for our purposes. 
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Figure 6. Displacements and fluxes for the case of Figure 5. Top left: contours of Re£ y , 
characteristic of the Alfven wave, and displaying phase mixing with increasing z along the 
resonance band. Top right: contours of F z Aif. Bottom left: F z f ast . Bottom right: shading 
representing absolute x-z wave energy flux \F XZ \ = (F x + F 2 ) 1 ' 2 . Different colour scaling 
ranges are used for the three panels, corresponding in each case to white being the maximum, 
and the broad background colour corresponding to zero. The maximum values for F z Aif, 
Fz, fasti and \F XZ \ are respectively 0.0044, 0.0072, and 0.0099. Fq = 1 is assumed throughout. 



exclusively plotting F Zi Aif- On the other hand, Fzfast (bottom left) quickly fades out 
in the z-direction. 

Finally, and most strikingly, the bottom right panel of Figure 6 displays contours 
of the total x-z flux, showing how it splits into the reflected and Alfvenic branches, 
each with lesser energy than the incident branch. One further interesting though sub- 
tle feature most apparent in this last graphic is that the outgoing ray has been shifted 
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Figure 7. The fluxes through x = — 2 (left panel) and the top boundary (right panel). In the 
latter case, F z Aif (solid) and F z f ast (dashed) arc displayed separately. 



slightly upwards relative to the ray theory path. We believe this is a consequence of 
the phase shift associated with the mode conversion. 

The reader is strongly advised to view the movie accompanying this paper as 
electronic supplementary material. It shows the dilatation field in greyscale (as 
in Figure 5a), but with the Alfven displacement £ y overplotted in blue and red, 
representing displacements out of and into the page respectively. The manner in 
which the two fields move in phase is most instructive. 

Figure 7 shows the energy balance between the incoming, reflected and Alfvenic 
ray bundles. Total incoming flux on x = —2, z < is 0.00761 and outgoing flux on 
x = -2, z > is 0.00417, with net flux 0.00343. The reflection coefficient is therefore 
0.55, which accords very well with 1 — Sk(a, 6) = 0.54 for the central wavenumber 
k z Q. The net outgoing flux on the top is 0.00340, matching almost perfectly (a small 
amount of energy is lost because of imperfect resolution on the grid and the low-order 
integration schemes used). 

By way of comparison, Figure 8 shows the total flux for the cases with k y = 1 and 
k y = 6 and the other parameters unchanged. These examples have central a = 0.053 
and 9 = 6.3° for k y = 1, and a = 1.505 and 6 = 33.7° for k y = 6, with corresponding 
narrower and wider evanescent gaps respectively and much lower central absorptions 
A = 0.12 and 0.063. As expected, the ray bundles suffer only very weak Alfven 
conversion (0.12 and 0.064), and there is little apparent upward displacement of the 
reflected rays. 

So how exactly are the reflection and absorption coefficients of a ray bundle related 
to those for the constituent modes? Integrating the physical-space x-flux F x {x,z) = 
F Re(ix*?a;) over all z at fixed x, it is easily shown that 



Ftot = / F x (x,z)dz = 



1 
2tt 



(Fx) 



2tt 



\A(k z )\ 2 F x (x,k z )dk z 



\A(k z )\ 2 dk z 



(22) 



(defining (. . .)), where F x is the flux in an individual mode k z , provided the in- 
tegrals exist (which is so for our gaussian amplitude). For the gaussian profile, 
2^ J_ \A(k z )\ 2 dk z = FIsJ^k. If we now split F x into its rightward and leftward 
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Figure 8. Same as the bottom left panel of Figure 6, except with k y = 1 (left) and k y 
(right). 



components at Xl, F tot = F x + F x , and similarly for F x , it is readily found that 



'K. 



<M 



F+ 



(23) 



where ^, is the overall reflection coefficient of the ray bundle and ^ is that for the 
individual Fourier modes. More simply, if _F+ is scaled to unity, then 1^ = (%). This 
shows how the reflection coefficient for a ray bundle may be constructed from those 
for the constituent modes. Clearly, the sharper is A in k z space, the more closely ^ 
matches ^. But this corresponds to a very broad ray bundle in physical space. If 
on the other hand we sharpen the bundle in z, it is broadened in k z , resulting on a 
broad flux- weighted average over k z . 

Selecting z well above the region where x is significant, Equation (21) may be 
integrated exactly to determine the Alfvenic flux at the top F Zitop (x) = F Zi au{x, oo) 
on each field line. Curiously, this physical-space quantity depends algebraically on 
the Fourier-space dilatation %: 



:,top 



■(*) 



^|A(fc(x))| 2 | X (x,fc(x))| 2 



(24) 



It is peaked close to where A(k(x)) is peaked, i.e., at k = u/a = fc z0 , the central 
resonant point, consistent with Figures 6 and 7. This expression for F z ^ op (x) is exact 
and does not suffer the slight numerical inaccuracies associated with the integrations 
for £. 
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Equation (22) and Equation (6) for F x may be combined to show 



Ftot(.x 2 ) - F tot (xi) = J F x (x 2 ,z)dz- / F x {x\,z)dz 

J -co J -co 

1 f'°° 
= — \A(k z )\ 2 (F x (x 2 ,k z )-F x (x 1 ,k z )) dk z (25) 

J — CO 

= F kl J" lA [ k ^f \X(X, fc(x))| 2 dx = J" F ZMJP dx , 

where x 2 > x\ and dk/dx = —k/2h has been used. Letting x 2 — > x^ confirms that 
indeed the energy entering and being consumed by the resonance at Xi is exactly 
that converted to F z along that field line. 

The significance of Equation (24) is that the Alfvenic flux is made up of precisely 
the modal dilatation on the resonances; at each point x the flux comes solely from 
the mode which has its resonance there. And providing A(k) is smooth, then so is 
F Zjt op, despite £ being undefined at the resonance. This is a pretty insight. 



6. Discussion and Conclusions 

The analysis of mode transmission and absorption presented here for the simple 
case of a cold MHD plasma allows us to at last understand the commonalities and 
distinctions between mode conversion and resonant absorption theories. 

Both the current model and the mode conversion models of Cally (2006) and 
Schunker and Cally (2006) share the characteristic that there is an avoided crossing 
or near passage of two distinct branches of the dispersion relation somewhere in x-k 
phase space, across which energy may tunnel. In the case of the Schunker and Cally, 
2006 analysis, an upward-travelling fast MHD wave with predominantly acoustic 
characteristics (since the sound speed c exceeds the Alfven speed a in the solar 
interior) is incident on a surface active region where a ~ c and transmits to a slow 
(but still predominantly acoustic) wave in the region above (a > c). The field-guided 
sound wave then takes the energy away to be either deposited in higher layers or to 
reflect from an acoustic cutoff. 

On the other hand, the modal resonant absorption model discussed here provides 
no opportunity for energy to be removed from the interaction region. Although a 
singular Alfven wave of sorts propagates along the z-directed field line at x = X, 
there is no net Alfvenic energy loss from any segment. Instead, energy is being 
continually stored deeper and deeper in the singularity (recall that it takes an infinite 
time for energy to actually reach X), or in reality is thermalized there by non-ideal 
effects. 

In particular, we may now understand resonant absorption as something that 
occurs when there is transmission to a "trapped" mode. This consideration deserves 
some further elaboration as it can easily be misinterpreted in various ways. 

First, one may have the impression that the resonances are related to the "closed- 
ness" of the system which artificially traps the energy at a certain position. This is 
strikingly untrue. In fact, in the modal model of Section 2 the field lines are open 
and hence an Alfvenic wave packet can propagate away in the z direction, though 
with no net Alfvenic loss, and yet a singularity still results. This brings us to the 



convRev.tex; 13/07/2010; 0:17; p. 15 



P.S. Cally and J. Andrics 



second possible misinterpretation, previously expressed by some authors (Schwartz 
and Bel, 1984; Hansen and Harrold, 1994), that the occurrence of resonances and the 
associated singularities is due to studying a too simple model (see Ruderman and 
Wright, 1999, for a slightly different slant on the Hansen and Harrold, 1994, results). 
A well-known example is the case of horizontal magnetic field in a gravitationally 
stratified atmosphere (e.g., Cally, 1984), where the linear wave equations reduce 
to a single second order ordinary differential equation with coefficient of the second 
derivative vanishing at Alfven and slow resonances thereby leading to classic resonant 
absorption. However, if the field is tilted only slightly, the equations become fourth 
order and the resonance is gone! The reason is that there has been mode conversion 
to a wave that can run upward along the field line to escape, whereas in the horizontal 
field case it has nowhere to go. This example, however, is misleading. The horizontal 
field case suffers the same shortcoming as our present model, i.e., the invariance 
along the field. The inclined field case seems to resolve that issue but at the cost of 
simplifying the model in another respect. In the case of an inclined field all field lines 
are alike and the propagation across the field is trivialised (by Fourier expansion). It 
is exactly this aspect that removes the resonances (Goossens, Poedts and Hermans, 
1985) which are essentially associated with the non-trivial behaviour of solutions 
across the field. Resonant absorption has indeed been studied in models including 
both variation along and across the field simultaneously (e.g., Thompson and Wright, 
1993). 

Based on these considerations it is clearly necessary to remove the invariance 
along field lines in order to progress further with our understanding. This is done 
in Section 5 by synthesizing a discrete ray bundle from the previously calculated 
singular normal modes. The resulting ray is not only smooth in \ as expected, 
but also smooth in the displacement £ as constructed directly from \. There is no 
sign in physical space of a local or even distributed energy build-up, and also the 
escaping Alfven wave is strikingly apparent. This Alfven wave excitation/conversion 
mediated by resonant absorption has been noted previously for kink waves in coronal 
loop models (Ruderman and Roberts, 2002). 

So is it "resonant absorption as mode conversion"? Or is it "mode conversion as 
resonant absorption"? It is both! We have fabricated a (mode converting) ray bundle 
from the (resonant absorbing) normal modes, but could just as well have constructed 
the normal modes from a collection of rays. 

The fundamental and inescapable conclusion therefore is that resonant absorption 
and mode conversion are indeed two aspects of the same process. The distinguishing 
characteristic is whether or not the system allows a net escape of energy from the 
absorbing region in physical space. If so, it is best understood in terms of mode con- 
version; if not, then energy escapes in wavenumber space only, producing arbitrarily 
small scales (or local non- ideal dissipation), and resonant absorption theory is the 
most natural description. 
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Appendix 

A. Numerical Solution for a General Profile 

A.l. Eikonal Description 

We follow the general asymptotic WKB expansion method as described by Bender 
and Orszag (1978), using e = fc" 1 as our "small parameter": 



X~exp 



fy- 1 ^ 



,n=0 



(26) 



and obtain a WKB hierarchy for So, 5*1, etc. Defining K 2 = f3 cos 2 6 — sin 2 = 

(P - Po) cos 2 e = n 2 - i, 

S' 2 = -K 2 (27a) 

2S' S[ + S{ - {p'/P)S' Q = (27b) 

S' 1 2 + 2S / Q S / 2 + S'{-(f3'/l3)S' 1 = (27c) 
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etc., whence 



S Q = 



Si = In 



y/-K{t) 2 dt 

(_ j8 )J(_ A: 2)-i 



(28a) 
(28b) 



Consequently, the physical optics approximation (i.e., up to and including S\) is 



X~(-f3)H-K 2 )-*{Aexp 



yJ-K{t) 2 dt 



+ -Bexp 



yJ-K(t) 2 dt 



(29) 
as e — > 0, where ^4 and B are integration constants. The validity of this approximation 
requires that eS^ <C 1 as e — > 0. In Region I, where both f3 and K 2 are positive, an 
alternate normalization will be used and So becomes imaginary, 



x~PHK 2 y 



Ai exp 



y/K(t)*dt 



+ B\ exp 



\/K(t) 2 dt 



(30) 



The two complex exponentials are respectively the incident and the reflected waves. 

A. 2. Frobenius Solution 

For the case of a general monotonic profile, Equation (3) is solved numerically using 
an initial value ODE solver starting from the decaying eikonal approximation (up to 
and including S2) at x 3> X, i.e., with A = and B = 1 in (29), and integrating left- 
ward. However, to avoid integrating through the singularity at x = X, the Frobenius 
solution (below) is used in a neighbourhood (X — S, X + 6). The numerical solution 
is then continued from X — S to a large negative x where the eikonal solution (29) 
may again be used to identify the incoming (A) and outgoing (B) parts. 

To develop the Frobenius solution about the resonance point assume that f3 is 
analytic in a neighbourhood of X, and write 



P = J2p n (x-X) n , 



(31) 



where /3 n = [^ n \X)/n\ and [i\ < 0. The dependent variable \ m ay then be expanded 
in a Frobenius series for which the indicial equation yields r = 2 and 0. The first 
solution is 



X 



a = ^ffln^- -X") 



n+2 



(32) 



whilst the second is 



Xb = X 



a \n(x-X) + ^A n (x-X) n . 



(33) 



The causal branch of the logarithm in this instance is the one where ln(a; — X) = 
ln(X — x) — m on x < X. 
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The coefficients are a = 1, and 



e" 2 (D p sin 2 6-C v cos 2 6) - E p 

( P 2 -l) ft < M > 



where 



p— 4 <j— 1 

C p = 2_^ a nlp-n-2 Iq = /, PkPq-k 
n=0 fc=l 

p-3 p-2 v 7 

Z3 P = 22 a n/3 P -,i-2 E p = 2_/n + 2)(2n + 1 - p)a n /3p - n 

n=0 n=0 

for p = 2, 3, .... The convention is adopted that a sum is zero if its upper limit is 
less than its lower limit. Thus for example D% = 0. 
Similarly, 

,2r> 1,2 

(36) 



(37) 













4> = - 


e 


2 Qi 


2e 2 

" sin 2 ' 






ft 


sin 2 (9 




Ai 


= A 2 


= 0. 


, and 




















Ap+l 


= 


(c|Vp T" {j~p 


+ e 


" 2 (^ 


sin 2 9 — L t 


, COS 2 0) 








(P 2 - 


1)01 




p = 


2,3, 


..., 


where 


















Vp 


p-i 

= £(*>- 

n=0 
p-1 


3n 


- 3)a„/3 p 


— n 


G p = 


n=l 
p-1 


2n + 3) 



(38) 

<Jp — / A-nPp—n Lip — / Si-n^fp— n • 

n=0 n=0 

With these formulae, the a n and A n coefficients may be built up recursively to any 
desired level. The leading terms of this expansion for the two linearly independent 
solutions are 

Xb = Xa ln(x - X) + 2, 2 esc 2 6 - (^ cot 2 + ||) (x - Xf + O ((* - X) 4 ) . 

(39) 
The logarithmic term is responsible for the absorption. 

A. 3. Reflection and Absorption Coefficients and Phase Retardation 

The overall effect of on an incident fast wave in Region I may be characterized by 
the reflection coefficient S^ = \Bi/Ai\ 2 and the associated absorption coefficient 
fA = 1 — %., as well as by the net phase retardation A<p = arg(y4.i_Bi). This may be 
achieved numerically for any chosen Alfven profile /?, but is primarily dependent on 
its local properties around the reflection/absorption region. 
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e r e * is the a; — > 
= e Sr e~ iSi that for 



It remains to justify the above expression for A(f>. If Xmc ~ < 
— oo asymptotics for the canonical incident wave and Xref ~ f 
the reflected wave, and \ = A Xinc + -S Xief is the combined solution which matches to 
exponential decay as x — > oo, then \ ~ eSr [ \A\e^ Si+ 't' A '> + \B\e~ i ( Si+( ~' t ' B ')'>~\, where 
4>A = arg A and <j>b = arg B. The incident wave is therefore advanced in phase by 4>a 
and the reflected wave is retarded by 4>b- The total phase retardation of the reflected 
wave with respect to the incident wave is therefore A0(er, 6) = 4>a + 4>b = arg AD. 



B. Linear Profile O 2 = 1 - x/h: The RAT Index 

Consider the case of the linear profile f2 2 = 1 — x/h. Letting £ = e~ 2 / 3 h~ 1 / 3 (x — X) 
and a = e -2 ' 3 /i 2 / 3 sin 2 0, Equation (3) takes the form (16). Subject to the boundary 
condition \ — > at £ — > oo and an arbitrary linear scaling, it has a unique solution 
X = 4>a for each value of the parameter a, and consequently a unique absorption 
coefficient, the RAT index, fiL(cr) which can be calculated numerically; it is tabulated 
in Table 1 as the 6 = column, displaying a peak of A a = 0.493698 at a = a = 
0.464434. Absorption is very well fitted by JZ(a) = 2.66904 <rexp[-0.363841cr 2 - 
1.81607ct], with maximum absolute error of 0.001 over ^ a ^ 3. 

A Frobenius expansion about the regular singular point £ = returns the two 
linearly independent solutions 



Xi 



71 = 



X2 = a Xl lnC + ^C„C" 

n=0 



(40) 



where negative £ should be identified with |£|e 17r . The coefficients are given by the 
recurrence formulae 



Co 



1, ci = 0, 



Cl 



Cre-3 + gC„_ 2 

n(n + 2) 



(41) 



and 



Co — 2, C\ — C2 — 0, C„ 



<rC n - 2 + Cn-3 ~ 2<?(n ~ l)c(n - 2) 
(n - 2)n 



(42) 



In the case a = the logarithmic term in X2 vanishes, leaving only Airy solutions 
Ai'(C) and Bi'(£), and indeed we expect the general solutions to resemble these in 
some sense for £ ^> a. 

As an alternative to the WKB solutions at large |x|, a local analysis about 
the irregular singular point £ = 00 may be used, revealing that the two linearly 
independent solutions ijj(a;C) and cj)(a;C,) have asymptotic behaviours 



C 1/4 exp 



± (f C 3/ 2 + aC l/2 



i^r /2 +- — 

4 s 32 



la A , o- 6 - 40o- 3 - 56 , „,, 

s 384 s 



+ 



(43) 
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as £ — > oo. We define ip(a; Q to be the "physical" solution for which tp — > as 
£ — > +00 along the positive real axis, i.e., 



^-C 1/4 exp 



f C 3/ 2 + aC i/2 



(44) 



for £ > 0. Note that, although they share the same contolling factor, ip(a; Q ^ f/)(0; £) 
as C — > oo, so a has a profound effect even at large £. For a = 0, (44) is just the 
asymptotic representation of — 2i/7r Ai'(£). 



C. Exponential Profile W 2 = eT x l h 



If Region III is sufficiently narrow (X small, i.e., 6 small) it suffices to replace j3 by its 
tangent there, but in practice this produces satisfactory accuracy only for \6\ < 15°. 
A better model is to replace h = — (dlogf2 2 /dx) _1 by its value at the resonance, 
i.e., to set fl 2 = e~ x / h with constant h. Then /? = e -^- x ^ h - 1, X = /ilnsec 2 6>, 
and P'(X) = -1/h. 

Table 1. Absorption coefficients Ji(cr,8) for the exponential profile. 



= 0° 



10° 



20" 



:-!0 J 



40° 



50° 



,.,! 



70° 



SO'-' 



90- 



0.0000 
0.0200 
0.0500 
0.1000 
0.2000 
0.3000 
0.4000 
0.4644 
0.5000 
0.6000 
0.7000 
. 8000 
0.9000 
1.0000 
1.1000 
1.2000 
1.3000 
1.4000 
1.6000 
1 . 8000 
2.0000 
2.2000 
2.4000 
2.6000 
2 . 8000 
3.0000 



0.0000 
0.0510 
0.1209 
0.2207 
0.3653 
0.4495 
. 4877 
0.4937 
0.4921 
0.4729 
0.4386 
0.3955 
0.3485 
0.3011 
0.2558 
0.2140 
0.1767 
0.1440 
0.0925 
0.0571 
0.0339 
0.0195 
0.0109 
0.0059 
0.0031 
0.0016 



0.0000 
0.0329 
0.1028 
. 2048 
. 3530 
. 4395 
0.4788 
. 4852 
. 4837 
. 4646 
. 4303 
0.3871 
0.3403 
0.2931 
0.2482 
0.2069 
0.1702 
0.1382 
0.0880 
. 0538 
0.0316 
0.0180 
0.0099 
0.0053 
0.0028 
0.0014 



0.0000 
0.0100 
0.0606 
0.1602 
0.3163 
0.4089 
0.4516 
0.4591 
0.4579 
0.4393 
0.4050 
0.3620 
0.3156 
0.2693 
0.2257 
0.1860 
0.1511 
0.1212 
0.0751 
0.0445 
0.0254 
0.0140 
0.0074 
0.0038 
0.0019 
0.0009 



0.0000 
. 0029 
0.0273 
0.1033 
0.2576 
0.3574 
0.4050 
0.4142 
0.4137 
0.3959 
0.3620 
0.3197 
0.2745 
0.2303 
0.1893 
0.1529 
0.1215 
0.0951 
0.0560 
0.0315 
0.0170 
0.0088 
. 0044 
0.0021 
0.0010 
0.0004 



0.0000 
0.0010 
0.0114 
0.0572 
0.1868 
0.2873 
. 3387 
. 3497 
. 3499 
. 3336 
0.3010 
0.2606 
0.2184 
0.1781 
0.1419 
0.1108 
. 0849 
0.0640 
. 0348 
0.0179 
. 0088 
0.0041 
0.0018 
. 0008 
. 0003 
0.0001 



0.0000 
0.0004 
. 0048 
. 0287 
0.1192 
0.2063 
0.2560 
0.2675 
0.2681 
0.2536 
0.2238 
0.1878 
0.1513 
0.1180 
0.0895 
0.0662 
. 0479 
. 0340 
0.0162 
0.0073 
0.0031 
0.0012 
0.0005 
0.0002 
0.0001 
0.0000 



0.0000 
0.0001 
0.0020 
0.0131 
0.0654 
0.1266 
0.1654 
0.1745 
0.1749 
0.1623 
0.1378 
0.1096 
0.0828 
0.0600 
0.0420 
0.0284 
0.0187 
0.0121 
. 0047 
0.0017 
0.0005 
0.0002 
0.0000 
0.0000 
0.0000 
0.0000 



0.0000 
0.0000 
0.0007 
0.0050 
0.0283 
0.0598 
0.0809 
0.0852 
. 0846 
. 0749 
0.0589 
. 0425 
. 0286 
0.0181 
0.0110 
. 0064 
0.0036 
0.0019 
0.0005 
0.0001 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 



0.0000 
0.0000 
0.0002 
0.0011 
0.0069 
0.0153 
0.0203 
0.0205 
0.0197 
0.0153 
0.0100 
0.0058 
0.0030 
0.0014 
0.0006 
0.0003 
0.0001 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 



0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
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Table 2. Phase retardations A<f>(a, 8) (degrees) for the exponential profile. The absence of 
several entries at a = 20 is due to numerical difficulties. 

a 6 = 0° 10° 20° 30° 40° 50° 60° 70° 80° 



0.0000 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.0200 


2.65 


2.58 


2.20 


1.79 


1.49 


1.29 


1.16 


1.07 


1.03 


0.0500 


6.71 


6.70 


6.44 


5.88 


5.29 


4.78 


4.41 


4.16 


4.02 


0.1000 


13.69 


13.74 


13.71 


13.37 


12.82 


12.21 


11.69 


11.30 


11.08 


0.2000 


28.36 


28.50 


28.80 


29.05 


29.12 


29.09 


29.04 


29.07 


29.20 


0.3000 


43.53 


43.74 


44.30 


45.03 


45.78 


46.54 


47.34 


48.22 


49.12 


0.4000 


58.45 


58.72 


59.47 


60.58 


61.90 


63.36 


64.97 


66.70 


68.38 


0.4644 


67.57 


67.87 


68.72 


70.01 


71.63 


73.47 


75.51 


77.70 


79.78 


0.5000 


72.36 


72.68 


73.58 


74.96 


76.71 


78.73 


80.98 


83.38 


85.63 


0.6000 


84.81 


85.15 


86.16 


87.73 


89.78 


92.19 


94.90 


97.75 


100.33 


0.7000 


95.65 


96.01 


97.08 


98.78 


101.02 


103.70 


106.70 


109.80 


112.46 


0.8000 


104.97 


105.35 


106.45 


108.23 


110.59 


113.42 


116.57 


119.74 


122.32 


0.9000 


112.97 


113.35 


114.47 


116.29 


118.71 


121.61 


124.79 


127.91 


130.30 


1.0000 


119.84 


120.23 


121.36 


123.18 


125.61 


128.51 


131.65 


134.62 


136.77 


1.1000 


125.79 


126.17 


127.30 


129.11 


131.51 


134.35 


137.37 


140.14 


142.04 


1.2000 


130.97 


131.34 


132.44 


134.22 


136.57 


139.32 


142.18 


144.70 


146 . 38 


1.3000 


135.49 


135.85 


136.93 


138.66 


140.93 


143 . 56 


146.23 


148.51 


149.98 


1 . 4000 


139.46 


139.82 


140.86 


142 . 53 


144.71 


147.18 


149.65 


151.71 


153.00 


1 . 6000 


146.07 


146.40 


147.36 


148.88 


150.83 


152.99 


155.07 


156.72 


157.75 


1 . 8000 


151.27 


151.56 


152.43 


153.80 


155.51 


157.36 


159.08 


160.43 


161.28 


2.0000 


155.40 


155.66 


156.44 


157.63 


159.12 


160.69 


162.13 


163.26 


163.97 


2.2000 


158.70 


158.94 


159.62 


160.67 


161.95 


163.29 


164.51 


165.47 


166.08 


2.4000 


161.37 


161.58 


162.18 


163.09 


164.20 


165.35 


166.40 


167.23 


167.77 


2.6000 


163.55 


163.73 


164.26 


165.05 


166.02 


167.02 


167.94 


168.67 


169.14 


2.8000 


165.34 


165.50 


165.96 


166.66 


167.51 


168.39 


169.21 


169.86 


170.28 


3.0000 


166.82 


166.97 


167.38 


168.00 


168.75 


169.54 


170.27 


170.85 


171.23 


20.0000 








179.31 


179.35 


179.39 


179.43 


179.47 


179.49 



Equation (15) is solved numerically (see Appendix A for details) subject to the 
boundary condition \ ~ ► as £ — > oo to find A\ and B\, and hence J? and A0, both 
regarded as functions of the independent parameters a and 0. These are plotted in 
Figure 2 and tabulated in Tables 1 and 2. 

The first panel in Figure 2 illustrates that A — > as a — > or oo at all angles 6. 
The peak occurs at a = a Q = 0.4644, A = 0.4937 for the RAT index 5\.{<j) = SL(a, 0) 
plotted here as a dashed curve. The peak at all other 6 is situated very close to a 
but reduces monotonically in height with increasing angle. 

But why this humped structure of A(a)? Recall that the distance X between 
the reflection and resonance points depends only on 6 (for a given profile fi 2 ), and 
not on e. As e decreases towards zero (a — > oo), the evanescent solution in Region 
III decays ever more quickly (see Equation 29), so energy tunnelling across the gap 
vanishes and with it %.. This is also consistent with our earlier observation that the 
fast and Alfven modes decouple as k y — > in Equation (1). On the other hand, as 
e — > oo (a — > 0), we again see a decay to zero. This is because the length scale of the 
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solution increases without bound and the Alfven singularity becomes progressively 
less important. This is seen explicitly in Equation (16) as a — > 0. 

The second panel of Figure 2 displays the phase retardation over the range ^ 
a ^ 3 and 0° ^ 6 ^ 85°, showing that it increases monotonically from at a = 
(k ± = 0) to 180° as a — > oo (k ± — > oo; this is confirmed by solution at much larger 
a) and depends only weakly on 9 (with a held fixed). 

C.l. ERAT Equation in Rational Form 

Equation (15) may be re-expressed as an equation with rational coefficients by 
making the further change of variable t = u(Q, resulting in 

(yt - \f X " +("-])*- K 1 " ^ + ff ] X = (45) 

which has regular singular points at t = and t = 1/v (x = X and x = +oo 
respectively) and an irregular singular point at t = oo. We have written cos 2 = 
1 — <jv here. Frobenius expansions about each of the regular singular points may be 
derived. 
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